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We perform a combination of searches for standard model Higgs boson production in pp colli- 
sions recorded by the DO detector at the Fermilab Tevatron Collider at a center of mass energy of 
^/s = 1.96 TeV. The different production and decay channels have been analyzed separately, with 
integrated luminosities of up to 9.7 fb~^ and for Higgs boson masses 90 < Mh < 200 GeV. We 
combine these final states to achieve optimal sensitivity to the production of the Higgs boson. We 
also interpret the combination in terms of models with a fourth generation of fermions, and models 
with suppressed Higgs boson couplings to fermions. The result excludes a standard model Higgs 
boson at 95% C.L. in the ranges 90 < Mh < 101 GeV and 157 < Mh < 178 GeV, with an expected 
exclusion of 155 < Mh < 175 GeV. In the range 120 < Mh < 145 GeV, the data exhibit an excess 
over the expected background of up to two standard deviations, consistent with the presence of a 
standard model Higgs boson of mass 125 GeV. 
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I. INTRODUCTION 

A fundamental goal of elementary particle physics is to 
understand the origin of clectroweak symmetry breaking. 
The proposed mechanism in the standard model (SM) in- 
troduces a doublet of complex scalar fields into the SM 
Lagrangian, the neutral component of which develops a 
vacuum expectation value that generates the longitudinal 
polarizations and masses of the W and Z bosons. This 
mechanism [1-4] gives rise to a single scalar boson, the 
Higgs boson (H), but does not provide a prediction for 
its mass. Fermions acquire their masses via their inter- 
actions with the scalar field. Precision clectroweak data, 
including the latest W boson and top quark mass mea- 
surements at the CDF and DO experiments at the Fermi- 
lab Tevatron Collider [5-7] , constrain the mass of a SM 
Higgs boson to Mh < 152 GeV [8] at 95% confidence 
level (C.L.). Direct searches at the ALEPH, DELPHI, 
L3, and OPAL experiments at the CERN e+e^ Collider 
(LEP) [9], the CDF and DO experiments [10, 11], and 
the ATLAS [12] and CMS [13] experiments at the CERN 
Large Hadron Collider (LHC) limit the SM Higgs boson 
mass to 122 GeV < Mh < 127 GeV at 95% C.L. The 
ATLAS and CMS Collaborations have each observed a 
new boson in its bosonic decay modes with a mass near 
125 GeV that is consistent with SM Higgs boson pro- 
duction [14, 15]. The CDF and DO CoUaborations have 
reported combined evidence for a particle consistent with 
the SM Higgs boson produced in association with a W 
or Z boson that decays to a bb pair [16]. 

In this Article, we combine the results of direct 
searches for the SM Higgs boson in pp collisions at y/s = 
1.96 TeV recorded by the DO experiment [17-20]. The 
analyses combined here search for signals of Higgs boson 
production through gluon-gluon fusion (GGF) {gg H) , 
in association with vector bosons {qq — >■ VH, where 
V = W,Z), and through virtual vector boson fusion 
(VBF) {qq — )• q'q'H). The analyses utilize data corre- 
sponding to integrated luminosities of up to 9.7 fb~^, 
collected during the years 2002-2011. The Higgs bo- 
son decay modes examined are H ^ bb, H ^ W~^W^, 
H — >■ r+r", and H — >■ 77. We organize the searches into 
analysis subchannels comprising different production, de- 
cay, and final state particle configurations, designed to 
maximize the sensitivity for each particular Higgs boson 
production and decay mode. 

We present an overview of the individual analyses in 
Section II. Section HI discusses the common methods of 
background estimation and simulation, while Section IV 
details the signal predictions and associated uncertainties 
used in the analyses. In Section V we describe the sta- 
tistical tehniques used in the combination, and provide 
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an overview of the most important systematic uncertain- 
ties. We validate our analysis techniques and statistical 
methods in Section VI by performing mesurements of the 
W Z + ZZ and WW production cross sections. In Sec- 
tion VII we present our results for the SM Higgs boson 
as well as two interpretations beyond the SM. We sum- 
marize our results in Section VIII. 



II. CONTRIBUTING ANALYSES 

A list of the analyses used in this combination is given 
in Table I. We summarize the analyses below, group- 
ing them according to the Higgs boson decay mode to 
which the analysis is most sensitive. To facilitate their 
combination, the analyses are constructed to be mutually 
exclusive after all event selections. 



A. H bb Analyses 

The most sensitive analyses for masses below Mh S 
130 GeV are those searching for H bb decays in associ- 
ation with a leptonically decaying V boson. To enhance 
the H ~> bb component in the data, the analyses use 
an algorithm (6-tagger) to identify jets that are consis- 
tent with 6-quark lifetime and fragmentation. Several 
kinematic variables sensitive to displaced vertices and to 
tracks with large transverse impact parameters relative 
to the production vertex are combined in a 6-tagging dis- 
criminant. This algorithm provides improvements when 
compared to the previously used artificial neural network 
(ANN) 6-tagger [31]. By adjusting the minimum require- 
ment on the output of the 6-tagger, a range of signal 
efficiencies and purities is achieved. 

The DO collaboration previously published a combina- 
tion oi H ^ bb analyses on the full Run II dataset [32]. 
The two searches focused on ZH production described 
below are unchanged from the previous combination, 
while the WH search differs slightly from the previous 
iteration in the multijet background estimation and a re- 
fined treatment of some systematic uncertainties. 

The WH Ivbb {£ = e, /i) analysis [21, 22] requires 
topologies with a charged lepton, significant imbalance 
in the transverse energy (-^y), and two or three jets (j). 
A boosted decision tree (BDT) discriminant [33-37] from 
TMVA [38] is used to discriminate against multijet back- 
ground. Using the average of the two highest outputs 
from the 6-tagger for all selected jets, six mutually ex- 
clusive 6-tagging categories are defined. Events with no 
6-tagged jets, and with exactly one of the lowest purity 
which can originate from a c quark in the hadronic decay 
W ^ cs are used for the H — >■ W~^W~ — > ii^q'q analysis, 
while the remaining events belong to the four 6-tagging 
categories that are used in the WH — )■ liybb analysis. A 
BDT discriminant is constructed for each lepton flavor. 
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TABLE L List of analysis channels, with the corresponding integrated luminosities and ranges in Mr considered in the combined 
analysis. See Section II for details. We group the analyses in four categories, corresponding to the Higgs boson decay mode to 
which the analysis is most sensitive: H ^ bh, H ^ W+W~ , H — > t'^t~ , and H — >■ 77. 



Channel {V ^W,Z and £ = e, ^i) 




Luminosity (fb ^) 


Mh (GeV) 


Reference 


WH — )> £ubb 




9.7 


90-150 


[21, 22] 


ZH ££bb 


H ^bb 


9.7 


90-150 


[23, 24] 


ZH — > wbb 




9.5 


100-150 


[25] 


H W+W- -> l+i^l-i? 




9.7 


100-200 


[26] 


H + X ^ W+W- fJ.^T^+ < 1 jet 




7.3 


155-200 


[27] 


H W+W- luq'q 
VH — eejj./ nne+X 
VH -> e+^l++X 


H W+W- 


9.7 
9.7 
9.7 


100-200 
100-200 
100-200 


[22] 
[28] 
[28] 


VH £pq'qq'q 




9.7 


100-200 


[22] 


VH ThThfl + X 

H + X^ £Thjj 


H T+r- 


8.6 
9.7 


100-150 
105-150 


[28] 
[29] 


— >■ 77 




9.7 


100-150 


[30] 



jet multiplicity, and 6-tagging category. In addition to 
kinematic variables, the inputs to the final discriminants 
include the 6-tagger output and the output from the mul- 
tijet discriminant. 

The ZH — )■ Ubb analysis [23, 24] requires two isolated 
charged leptons and at least two jets, at least one of 
which must pass a tight ^-tagging requirement. A kine- 
matic fit corrects the measured jet energies to their best 
fit values according to the constraints that the dilepton 
invariant mass should be consistent with the Z boson 
mass Mz and the total transverse momentum of the lep- 
tons and jets should be consistent with zero. The events 
are divided into "double-tag" and "single-tag" subchan- 
nels depending on whether a second jet passes a loose 
6-tagging requirement. The analysis uses random forest 
(RF) [38] discriminants to provide distributions for the 
final statistical analysis, applied in a two-step process. 
First, the events are divided into independent it-depleted 
and tt-enriched subchannels using a dedicated RF that 
is trained to discriminate signal from the tt backgrounds 
in each lepton and 6-tagging subchannel. Final discrim- 
inants are then constructed to separate signal from all 
backgrounds. The limit is calculated using the output 
distributions of the final discriminants for both the tt- 
depleted and tt-enrichcd samples. The H + X iT^jj 
analysis, where denotes r-lepton decays into hadrons, 
discussed in Sec. II C includes a contribution from ZH 
production with Z — r+r" and H ^ bb decays. 

The ZH — > i^Dbb analysis [25] selects events with large 
and two jets. This search is also sensitive to the WH 
process when the charged lepton from W (ly decay 
is not identified. Events selected in the WH — )• £vbb 
analysis are rejected to ensure no overlap between the 
two analyses. About 47% of signal in this analysis comes 
from WH — > ii/bb events in which the charged lepton 
fails the WH — iiybb analysis selection requirements. 
Variables such as significance and a track-based miss- 
ing transverse momentum are used to reject events with 



]^rp arising from mismeasurement of jet energies. The 
multijet background is further reduced by employing a 
dedicated BDT discriminant before applying 6-tagging. 
Two 6-tagging subchannels are defined using the sum of 
the 6-tagging discriminant outputs of the two jets. BDT 
classifiers, trained separately for different 6-tagging cat- 
egories, are used as a final discriminant. 



B. H ^ WW* Analyses 

We search for Higgs boson decays to two W bosons 
from the three dominant production mechanisms: gluon- 
gluon fusion, associated production, and vector boson 
fusion. In H W+W' decays with Mh < 2Mw, at 
least one of the W bosons will be virtual {W*). 

The dominant search channels are H — > W+W' — ?> 
{e+ e- , fj.+ fj- , IJ,^)^!? [26]. The presence of neutrinos 
in the final state prevents precise reconstruction of the 
candidate Mh- Events are characterized by large I^t and 
two isolated leptons of opposite electric charge. Each fi- 
nal state is further subdivided according to the number 
of jets in the event: no jets, one, and more than one jet. 
This division requires an evaluation of theoretical uncer- 
tainties on the signal predictions for each jet category, as 
will be discussed in Section IV. 

The dielectron and dimuon channels use BDT discrim- 
inants to reduce the dominant Drell-Yan background, 
while the e^/.i^ channel uses ^j^-related variables to min- 
imize backgrounds. All channels separate events into 
WM^-enriched and M^VF-depleted subchannels. In the 
dielectron and dimuon channels, dedicated BDTs are ap- 
plied to events with no jets or exactly one jet. Events 
with no jets are split according to the lepton quality in 
the e^fi+ channel. BDT response distributions, using 
several kinematic variables as inputs, are used as final dis- 
criminants. Inputs also include 6-tagging information for 
subchannels containing jets to reject the tt background. 
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We consider final states where at least one W boson 
decays to tv, and the r lepton decays into hadrons (t/j) 
and Vr [H + X ^ W+W' -> /^=^t^+ < 1 jet) [27]. 
Final states involving other t decays and misidentified 
Th decays are included in the H — > W^W~ analyses 
channels. This channel uses ANN outputs [38] for a final 
discriminant. 

We also include analyses that search for H — > W^W~ 
with one of the W bosons decaying into qq'. The H ^• 
W^W~ — lv(lq analysis [22] has the same initial selec- 
tions as the WR — ?> Ivhb search, except that it considers 
only events with no 6-tagged jets, and with exactly one 
6-tagged jet of the lowest purity that can originate from 
a c quark. The RF discriminants trained for each lepton 
fiavor, jet multiplicity, and 6-tagging category serve as 
the final discriminant variables. 

For VH — > VWW production, we consider final 
states containing: (i) three charged Icptons {VH ^• 
eefj,/ fifie+X) [28]; (ii) an electron and muon with the 
same charge (e^/i^ + X) [28]; and (iii) final states 
with one lepton, and at least four jets {VH — > 
iuq'qq'q) [22]. 

The VH — )- eefi/ fifie+X analyses use BDT outputs as 
final discriminants. In the ^/xe final state, events are split 
into three mutually exclusive regions to separate signal 
from Z+jets and other backgrounds. 

The e^/i^ + X analysis, in which the same-sign re- 
quirement suppresses the Drell-Yan background, uses a 
two-step multivariate approach: (i) a BDT is used to sup- 
press most of the dominant backgrounds from multijet, 
Vl^-l-jets, and W + j events, and (ii) another BDT is used 
to discriminate signal from the remaining backgrounds. 

The VH — >■ Ivq'qq'q analysis [22] has selections similar 
to the H — >■ W~^W~ — >■ ivq'q analysis, but requires at 
least four jets. Separate BDTs are trained for different 
backgrounds, and then they are used as input variables 
to the final RF discriminant. 



C. H ^ T+T- and H 77 analyses 

Higgs boson decays involving r leptons are included in 
different ways. The VH — > T^ThH + X analysis [28] uses 
a two-stage BDT approach, in which the first BDT dis- 
criminates between signal and backgrounds other than 
diboson (VV) production, and the second BDT, trained 
to distinguish between signal and all backgrounds, is im- 
plemented after selecting events that pass the first BDT 
requirement. 

The H + X (-Thjj analysis [29] selects events with 
one electron or muon, a r/j, and two or more jets. It 
is sensitive to associated VH, VBF, and gg ^ H + X 
production, and to both H ^ tt and H WW de- 
cays. A BDT, trained to distinguish between signal with 
H ^ TT and H — >• WW decays, is used to create rr- 
and W^W^-dominated subchannels within the electron and 



muon channels. Each of the four resulting subchannels 
has a BDT as the final discriminant. 

We also include in the combination an analysis that 
searches for Higgs boson decaying to two photons [30]. 
The Higgs boson is assumed to be produced via GGF, 
VBF, and associated VH production. The contribution 
of jets misidentified as photons is reduced by combining 
information sensitive to differences in the energy deposi- 
tion in the tracker, calorimeter and central preshower in 
an ANN for each photon candidate. The ANN output de- 
fines photon-dominated and jet-dominated regions, each 
of which is split into signal-rich and signal-depleted con- 
tributions based on the diphoton invariant mass. A BDT 
built with ten variables, including the diphoton mass, 
serves as the final discriminant in the signal rich region, 
while the diphoton mass only is the final discriminant in 
signal-depleted regions. 



III. BACKGROUND ESTIMATION 

All analyses estimate backgrounds from multijet pro- 
duction through special data control samples. The 
other backgrounds are determined from Monte Carlo 
(MC) simulation. MC samples are generated us- 
ing the PYTHIA [39], ALPGEN [40], SHERPA [41], Or 

SINGLETOP [42, 43] event generators, with pythia 
also providing parton showering and hadronization for 
ALPGEN and SINGLETOP. All generators use the 
CTEQ6L1 [44, 45] leading order (LO) parton distri- 
bution functions (PDF). Drell-Yan and PF-|-jets yields 
are normalized to next-to-next-to-LO (NNLO) calcula- 
tions [46], or, in some analyses, to data control sam- 
ples [23, 24, 26, 28]. For the V + bb/cc MC samples, 
generated separately from the l^-|-light-flavor events, we 
apply additional normalization factors calculated at next- 
to-LO (NLO) from mcfm [47, 48] to account for the 
heavy-flavor to light-flavor production ratio. Diboson 
background cross sections are normalized to NLO calcu- 
lations from MCFM. Top quark pair and single top quark 
production are normalized to approximate NNLO [49] 
and next-to-NNLO (NNNLO) [50] calculations, respec- 
tively. We correct the transverse momentum (pt) spec- 
trum of the Z boson in the MC to match that observed in 
data [51]. We correct the W boson using the same de- 
pendence, taking into account differences between the pt 
spectra of Z and W bosons predicted in NNLO QCD [52] . 
We account for Wj* production and its interference with 
WZ production using powheg [53] in analyses where 
this effect is significant: H W^W~ evev j [W[iv , 
VH e^p.'^+X, and VH ee^i/ ^i^ie+X . 
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IV. SIGNAL PREDICTIONS AND 
UNCERTAINTIES 

An outline of the procedures for the signal predic- 
tions and associated uncertainties is given below. Ref- 
erence [10] contains a more complete discussion. 

We simulate signal with pythia using the CTEQ6L1 
PDFs to model the parton shower, fragmentation, and 
hadronization. We reweight the Higgs boson px spec- 
tra for GGF production to the prediction obtained from 
HQT [54-56]. To evaluate the impact of the scale un- 
certainty on the differential spectra, we use the RES- 
BOS [57, 58] generator and apply the scale-dependent 
differences in the Higgs boson px spectrum to the HQT 
prediction. We propagate these changes to the final dis- 
criminants as a systematic uncertainty on the differential 
distribution which is included in the calculation of the 
limits. 

We normalize the Higgs boson signal predictions to 
the most recent higher-order calculations (see Table H). 
The gg — > H production cross section (agg^n) is calcu- 
lated at NNLO in QCD with a next-to-next-to-leading- 
log resummation of soft gluons. The calculation also 
includes two-loop electroweak effects and the running 
&-quark mass [59, 60]. The values in Table H arc up- 
dates [61] of these predictions, with the top quark mass 
set to 173.1 GeV [62], and includes an exact treatment 
of the massive top quark and bottom quark loop correc- 
tions up to NLO and next-to-leading-log (NLL) accuracy. 
The factorization scale and renormalization scale fXR 
choices for this calculation are hf = A*fl = ^h- These 
calculations are improvements over the previous NNLO 
calculations of (jgg^H [63-65] . We apply the electroweak 
corrections computed in Refs. [66, 67]. The soft gluon re- 
summation uses the calculations of Ref. [68] . The gluon 
PDF and the accompanying value of as{q^) strongly in- 
fluence Ugg^H- The cross sections we use are calculated 
with the MSTW 2008 NNLO PDFs [69], as recommended 
by the PDF4LHC working group [70, 71]. 

For analyses that consider inclusive gg ^ H produc- 
tion, but do not split the signal into separate channels 
based on the number of reconstructed jets, we use the 
uncertainties on inclusive production from the simulta- 
neous variation of the factorization and renormalization 
scale up and down by a factor of two. We use the pre- 
scription of the PDF4LHC working group for evaluating 
PDF uncertainties on the inclusive production cross sec- 
tion. QCD scale uncertainties that affect the cross section 
via their impact on the PDFs are included as a correlated 
part of the total scale uncertainty. The remainder of the 
PDF uncertainty is treated as uncorrclated with the un- 
certainty on the QCD scale. 

For analyses of gg ^ H production that divide events 
into separate channels based on the number of recon- 
structed jets, we evaluate the impact of the scale uncer- 
tainties following the procedure of Ref. [72]. We treat 
as uncorrelated the QCD scale uncertainties obtained 



from the NNLL inclusive [59, 60], NLO with one or more 
jets [73], and NLO with two or more jets [74] cross section 
calculations. We then obtain QCD scale uncertainties for 
the exclusive gg — >■ 7J-|-n jets (n = 0, 1, > 2) categories by 
propagating the uncertainties on the inclusive cross sec- 
tion predictions through the subtractions needed for the 
exclusive rates. For example, wc obtain the H+0 jet cross 
section by subtracting the NLO H+ > 1 jets cross sec- 
tion from the inclusive NNLL-f-NNLO cross section. We 
therefore assign three separate, uncorrelated QCD scale 
uncertainties that lead to correlated and anticorrelated 
contributions between exclusive jet categories. The pro- 
cedure in Ref. [73] is used to determine the uncertainties 
from the choice of PDF. These are obtained separately 
for each jet bin and treated as fully correlated between 
jet bins. 

Another source of uncertainty in the prediction of 
Ugg^H is the extrapolation of QCD corrections computed 
for heavy top-quark loops to the light-quark loops in- 
cluded as part of the electroweak corrections. Uncertain- 
ties at the level of 1-2% are already included in the cross 
section values we use [59, 60]. The factorization of QCD 
corrections is expected to be reliable for Mh values much 
larger than the masses of the particles contributing to the 
loop [59]. There is a 4% change in the predicted cross 
section when removing all QCD corrections from the di- 
agrams containing light-flavored quark loops. For the 
6-quark loop [59], the QCD corrections are much smaller 
than for top-quark loops, confirming that the procedure 
does not introduce significant uncertainties. We there- 
fore do not consider any additional uncertainties from 
this source. 

For WH and ZH production we use cross sections 
computed at NNLO [75]. This calculation starts with 
the NLO calculation of v2hv [76] and includes NNLO 
QCD contributions [77], as well as one-loop electroweak 
corrections [78]. For VBF production, we use the VBF 
cross section computed at NNLO in QCD [79]. Elec- 
troweak corrections to the VBF production cross section, 
computed with the hawk program [78] are included al- 
though they are very small (< 0.03 fb) for the Mh range 
that wc consider. 

The predictions of Higgs boson decay branching frac- 
tions, B, are taken from hdecay [80, 81], and are also 
listed in Table H. Uncertainties on the branching frac- 
tions are taken from Ref. [82] . 



V. LIMIT CALCULATIONS 

We combine results using the CLs method with a neg- 
ative log-likelihood ratio (LLR) test statistic [83, 84] for 
the signal-plus-background (s -I- h) and background-only 
{h) hypotheses, where LLR = — 21n(Ls+f,/Lf,), and Lhy is 
the likelihood function for the hypothesis hy. The value 
of CLs is defined as CLg = CLg+b/CLi, where CL^+b 
and CLb are the confidence levels for the s + b and the 
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TABLE II: The production cross sections (in fb) and decay branching fractions (in %) for each SM Higgs boson mass considered 
in tire combined analysis. 



Mh (GeV) 






0"ZH 


0"VSF 


B{H -> bb) 


B{H cc) 


B{H T+T-) 


B{H W+W-) 


B(H ZZ) 


B{H 77) 


90 


2442 


394.7 


224.0 


118.2 


81.2 


3.78 


8.41 


0.21 


0.042 


0.123 


95 


2101 


332.1 


190.3 


108.8 


80.4 


3.73 


8.41 


0.47 


0.067 


0.140 


100 


1821 


281.1 


162.7 


100.2 


79.1 


3.68 


8.36 


1.11 


0.113 


0.159 


105 


1584 


238.7 


139.5 


92.3 


77.3 


3.59 


8.25 


2.43 


0.215 


0.178 


110 


1385 


203.7 


120.2 


85.2 


74.5 


3.46 


8.03 


4.82 


0.439 


0.197 


115 


1215 


174.5 


103.9 


78.7 


70.5 


3.27 


7.65 


8.67 


0.873 


0.213 


120 


1072 


150.1 


90.2 


72.7 


64.9 


3.01 


7.11 


14.3 


1.60 


0.225 


125 


949 


129.5 


78.5 


67.1 


57.8 


2.68 


6.37 


21.6 


2.67 


0.230 


130 


842 


112.0 


68.5 


62.1 


49.4 


2.29 


5.49 


30.5 


4.02 


0.226 


135 


750 


97.2 


60.0 


57.5 


40.4 


1.87 


4.52 


40.3 


5.51 


0.214 


140 


670 


84.6 


52.7 


53.2 


31.4 


1.46 


3.54 


50.4 


6.92 


0.194 


145 


600 


73.7 


46.3 


49.4 


23.1 


1.07 


2.62 


60.3 


7.96 


0.168 


150 


539 


64.4 


40.8 


45.8 


15.7 


0.725 


1.79 


69.9 


8.28 


0.137 


155 


484 


56.2 


35.9 


42.4 


9.18 


0.425 


1.06 


79.6 


7.36 


0.100 


160 


432 


48.5 


31.4 


39.4 


3.44 


0.159 


0.397 


90.9 


4.16 


0.0533 


165 


383 


43.6 


28.4 


36.6 


1.19 


0.0549 


0.138 


96.0 


2.22 


0.0230 


170 


344 


38.5 


25.3 


34.0 


0.787 


0.0364 


0.0920 


96.5 


2.36 


0.0158 


175 


309 


34.0 


22.5 


31.6 


0.612 


0.0283 


0.0719 


95.8 


3.23 


0.0123 


180 


279 


30.1 


20.0 


29.4 


0.497 


0.0230 


0.0587 


93.2 


6.02 


0.0102 


185 


252 


26.9 


17.9 


27.3 


0.385 


0.0178 


0.0457 


84.4 


15.0 


0.00809 


190 


228 


24.0 


16.1 


25.4 


0.315 


0.0146 


0.0376 


78.6 


20.9 


0.00674 


195 


207 


21.4 


14.4 


23.7 


0.270 


0.0125 


0.0324 


75.7 


23.9 


0.00589 


200 


189 


19.1 


13.0 


22.0 


0.238 


0.0110 


0.0287 


74.1 


25.6 


0.00526 



b hypotheses, respectively. These confidence levels are 
evaluated by integrating the corresponding LLR distribu- 
tions populated by simulating outcomes assuming Pois- 
son statistics. Separate channels and bins are combined 
by summing LLR values over all bins and channels. This 
method provides a robust means of combining channels 
while maintaining each individual channel's sensitivity 
and different systematic uncertainties. Systematic un- 
certainties are treated as nuisance parameters with Gaus- 
sian probability distributions constrained by their priors. 
This approach ensures that the uncertainties and their 
correlations are propagated to the outcome with their 
appropriate weights. 

To minimize the degrading effects of systematic uncer- 
tainties on the search sensitivity, we fit the individual 
background contributions to the observed data by maxi- 
mizing a likelihood function [85] . The likelihood is a joint 
Poisson probability over the number of bins in the calcu- 
lation and is a function of the nuisance parameters and 
their uncertainties. The maximization of the likelihood 
function is performed over the nuisance parameters, with 
separate fits performed to both the b and s + b hypotheses 
for each Poisson MC trial. We have verified that all fit 
parameters and pulls on the systematic uncertainties are 
well-behaved. 

The CLg approach used in this combination utilizes 
binned final variable distributions rather than a single 
fully integrated value for each contributing analysis. The 
signal exclusion criteria are determined by increasing the 
signal cross section until CLg < 0.05, which defines a 
signal cross section excluded at the 95% C.L. 



A. Final Variable Distributions 

All analyses are performed for the Mh range listed 
in Table I at 5 GeV intervals. Each analysis provides 
binned distributions of its final discriminants for each 
value of Mh and subchannel. The input distributions 
for individual channels can be found in the corresponding 
references in Table I. 

The limit calculation uses the full information avail- 
able in the individual discriminants. However, for visu- 
alization purposes it can be useful to collect all of the 
inputs into a single distribution. To preserve sensitivity 
from the bins with high signal-to-background {s/b) ra- 
tios, where s is the number of signal and b the number 
of background events, only bins with similar s/b ratio 
are combined. The aggregate distribution is formed by 
reordering all of the bins from the input distributions ac- 
cording to s/b ratio. The range of s/b ratio is large, so 
logiQ{s/b) is used. Figure 1 shows the aggregate distri- 
butions for Mh = 125 GeV and Mh = 165 GeV, in- 
dicating good agreement between data and predictions 
over several orders of magnitude. Figure 2 shows the 
same distributions after subtracting the expected back- 
ground from the data, where solid lines represent the ±1 
standard deviations (s.d.) in systematic uncertainty af- 
ter a fit to the background-only hypothesis. Integrating 
the distributions in Fig. 1 from the highest to the lowest 
s /b events illustrates how the data compare to the b and 
s + b hypotheses as the events in the highest s/b bins ac- 
cumulate. Figure 3 shows these cumulative distributions 
for approximately 150 of the most significant events as a 
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function of the integrated number of signal events. For 
Mh = 125 GeV, the highest s/h bins contain an excess 
of signal-hke events, while for Mh — 165 GcV, the data 
follow the background-only expectation. 



B. Systematic Uncertainties 

Systematic uncertainties on signal and backgrounds 
vary among the analyses and they are described in detail 
in Refs. [21-30]. We summarize below only the major 
components. Most analyses have an uncertainty of 6.1% 
from the integrated luminosity [86] , while the overall nor- 
mahzations in the ZH IM, H W+W' i+iyi'D, 
and VH e^i^i^+X analyses are determined from the 
mass peak oi Z ^ ££ and Z — >■ r+r" decays in data 
assuming the NNLO Zf-f* cross section, reducing the 
uncertainty to about 1%. The H ^ bb analyses have 
an uncertainty of 1-10% due to the uncertainty on the 
6-tagging rate, depending on the number and quality of 
tagged jets. All analyses take into account uncertainties 
on jet-energy scale, resolution, and jet identification effi- 
ciency, for a combined uncertainty of w 7%. All analyses 
include uncertainties associated with measurement and 
acceptances of Icptons, which range from 1% to 9% de- 
pending on the final state. The largest contribution to all 
analyses is from the uncertainty on the simulated back- 
ground cross sections which are 4-30% depending on the 
specific background process. These values include both 
the uncertainty on the theoretical cross section calcula- 
tions and the uncertainties on the higher-order correc- 
tion factors. The uncertainty on the expected multijet 
background in each channel is dominated by the statis- 
tics of the data sample from which it is estimated. It is 
considered separately from the uncertainties on the sim- 
ulated backgrounds' cross sections, and ranges from 10% 
to 30%. All analyses take into account the uncertainties 
on the differential cross sections arising from the choice of 
PDF set and QCD scale. The H ^ W+W' ^ l+ul'v 
{£ — e, fi) analyses divide the data according to jet mul- 
tiplicity, and consider uncertainties on the contribution 
from GGF that are a function of jet multiplicity. In addi- 
tion, several analyses incorporate uncertainties that alter 
differential distributions and kinematics of the dominant 
backgrounds in the analyses. These uncertainties are es- 
timated from the variation of the final discriminant dis- 
tribution due to generator and background modeling un- 
certainties. Correlations between systematic sources are 
also carried through in the calculations. For example, 
the uncertainty on the integrated luminosity is taken to 
be fully correlated between all signals and backgrounds 
obtained from simulation. Hence any fluctuation in lu- 
minosity is common to all channels for a single pseudoex- 
periment. All systematic uncertainties originating from 
a common source are assumed to be fully correlated. 



VI. ANALYSIS TECHNIQUE VALIDATION 
WITH DIBOSON PRODUCTION 

To validate our analyses techniques, we measure di- 
boson production cross sections in the V + bb and £i^£iy 
final states. The analyses use multivariate discriminants 
that utilize the same input variables as the discriminants 
used for the Higgs boson search, but with one or more 
diboson processes acting as the signal. The modified 
WH £M, ZH ££bb, and ZH uvbb analyses 
(collectively called the V Z analyses) treat the W Z and 
ZZ processes as signal, and the WW process as a back- 
ground. The Higgs boson processes are not taken into 
account in this validation procedure. Figure 4(a) shows 
the background-subtracted data for the dijet invariant 
mass in the VZ analyses, and Fig. 4(b) for the combined 
output of the V Z discriminant. Similarly, the modified 
H — > W^W~ — > £^i'£~D analysis uses the WW process 
as the signal with the WZ and ZZ processes as back- 
grounds. Figure 4(c) shows the background-subtracted 
data for the output of the WW discriminant. The VZ 
analyses measure a WZ + ZZ production cross section 
of 0.73 ±0.32 times the SM prediction of 4.4 pb obtained 
with MCFM. The significance for this measurement to be 
non-zero is 2.4 s.d. with an expected significance of 3.4 
s.d. The WW production cross section is measured to 
be 1.01 ± 0.06 times the SM prediction of 11.3 pb, also 
based on MCFM. Both measurements confirm our ability 
to extract a small signal from a large background in the 
same final states, using the same analysis techniques as 
the search for the Higgs boson, providing validation of 
the background modeling. 



VII. HIGGS BOSON RESULTS 

A. Limits on standard model Higgs boson 
production 

We obtain limits on the product of the Higgs boson 
production cross section, an, and branching fractions 
B{H — > 66/VF+W^~/t+t^/77) using individual chan- 
nels [21-30]. We present results in terms of the ratio 
of the upper limit on (Th at 95% C.L. relative to the SM 
predicted values as a function of Mh, where the rela- 
tive cross sections and branching fractions are kept as 
predicted by the SM. The SM prediction is therefore ex- 
cluded at the 95% C.L. for the Mh values at which the 
ratio falls below unity. 

The LLR distributions for the full combination are 
shown in Fig. 5. Included in these figures are the median 
LLR values expected for the s + b hypothesis (LLR^+f,), 
b hypothesis (LLRb), and the results observed in data 
(LLRobs)- The shaded bands represent the ±1 and 
±2 s.d. departures for LLR;,. These distributions can 
be interpreted as follows: 
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FIG. 1: (color online) Distributions of logj^o(s/&) for data from all contributing channels for (a) Mr = 125 GeV, and (b) 
Mh ~ 165 GeV after a fit to data assuming the background-only hypothesis. The data (points with Poisson statistical errors 
on the expected number of signal+background events) are compared to the expectation from background (light shaded) and 
signal (dark shaded). 
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FIG. 2: (color online) Background-subtracted distributions as a function of logj^Q(s/6) for data from all contributing channels 
for (a) Mh = 125 GeV and (b) AIh ~ 165 GeV after a fit to data assuming the background-only hypothesis. The background- 
subtracted data (points with Poisson statistical errors on the expected number of signal-|-background events) are compared to 
the expected signal (shaded). The solid lines represent the ±1 s.d. systematic uncertainty on the background after the fit. 
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FIG. 3: (color online) Cumulative number of events after integrating the final discriminant bins in decreasing order of logjQ(s/b) 
until the expected signal yield indicated on the x-axis is reached for (a) Mr = 125 GeV and (b) Mr = 165 GeV. The integrated 
b-only and s + b predictions are shown after their respective fits as a function of the accumulated number of signal events. The 
points show the integrated number of observed events with statistical errors. 
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FIG. 4: (color online) Background-subtracted data (points 
with statistical errors), measured diboson signal, and system- 
atic uncertainties after a fit to the s + b hypothesis for (a) 
the dijet invariant mass in the combined VZ — >■ Vbb analy- 
ses, (b) the output of the multivariate discriminant for the 
VZ — >■ Vbb analyses, rebinned in logjQ s/b, and (c) the out- 
put of the multivariate discriminant for the WW analysis, 
rebinned in log^Q s/b. The solid lines represent the ±1 s.d. sys- 
tematic uncertainty constrained by data. 



(i) The separation between LLRf, and LLRs_|_6 provides 
a measure of the discriminating power of the search, 
and illustrates the effectiveness of the analysis to 
separate the s + b and b hypotheses. 

(ii) The width of the LLRj, distribution (shown here 
as ±1 and ±2 s.d. bands) provides an estimate of 
the sensitivity of the analysis to a signal-like back- 
ground fluctuation in the data, taking the system- 
atic uncertainties into account. For example, the 
sensitivity is limited when a 1 s.d. background fluc- 
tuation is large compared to the difference between 
the s + b and b expectations. 

(iii) The value of LLRobs relative to LiLRs+b and LLRf, 
indicates whether the data distribution appears to 
be more s + &-likc or 6-like. The significance of any 
departures of LLRobs from LLR;, can be evaluated 
through the width of the LLRf, distribution. 

As shown in Table I, only the WH — >■ ii/bb and 
ZH — 7> ££bb channels contribute to the combination below 
Mh = 100 GeV. Figure 5 shows that the observed LLR is 
compatible with the s + b hypothesisfor 120 < Mh < 145 
GeV. 

Figure 6 shows the expected and observed upper lim- 
its on an at 95% C.L. relative to the SM, for the mass 
region 90 < Mh < 200 GeV, for all analyses combined. 
These results arc also summarized in Table III. We ex- 
clude the SM Higgs boson at 95% C.L. in the mass ranges 
90 < Mh < 101 GcV and 157 < Mh < 178 GeV. Our 
expected exclusion range is 155 < Mh < 175 GcV. 
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FIG. 5: (color online) The observed (black solid line) and 
expected LLRs for the b (black short-dashed line) and s -f & 
hypotheses (red/light short-dashed line), as well as the LLR 
expected in the presence of a SM Higgs boson with Mh ~ 
125 GeV (blue long-dashed line) for all analyses combined 
for the range 90 < Mh < 200 GeV. The shaded bands 
correspond, respectively, to the regions enclosing ±1 and 
±2 s.d. fluctuations of the background. 



Figure 7 shows the values for the observed CLs+b and 
its expected behavior as a function of Mh ■ The quantity 
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TABLE III: Expected (median) and observed upper limits on the cross sections relative to the SM at 95% C.L. for the combined 
analyses for the range 90 < Mh < 200 GeV. 





90 


95 


100 


105 


110 


115 


120 


125 


130 


135 


140 


145 


150 


155 


160 


165 


170 


175 


180 


185 


190 


195 


200 


Expected 
Observed 


1.29 
0.96 


1.40 
0.89 


1.13 
0.95 


1.21 
1.39 


1.32 
1.39 


1.45 
1.99 


1.59 
2.66 


1.66 
2.92 


1.69 
2.56 


1.58 
2.79 


1.49 
2.88 


1.33 
2.36 


1.17 
1.84 


1.02 
1.23 


0.75 
0.78 


0.70 
0.66 


0.86 
0.75 


1.02 
0.85 


1.21 
1.11 


1.55 
1.31 


1.89 
1.96 


2.22 
2.85 


2.55 
3.12 
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FIG. 6: (color online) Expected (median) and observed ratios 
for the upper limits of the cross section an at 95% C.L. rela- 
tive to the SM values for all analyses combined for the range 
90 < Mh < 200 GeV. The shaded bands correspond to the re- 
gions enclosing ±1 and ±2 s.d. fluctuations of the background, 
respectively. The long-dashed line represents the expectation 
if a Mh ~ 125 GeV Higgs boson were present in the data 
with the SM cross section. 



FIG. 7: (color online) The observed (black solid line) and 
expected CLs-i-b {s + b p-value) for the no-Higgs boson hy- 
pothesis (black short-dashed line) and in the presence of a 
SM Higgs boson with Mh = 125 GeV (blue long-dashed line) 
for all analyses combined for the range 90 < Mh < 200 GeV. 
The shaded bands correspond, respectively, to the regions en- 
closing ±1 and ±2 s.d. fluctuations of the background. The 
three red horizontal dashed lines indicate the p-values corre- 
sponding to significances of 1, 2 and 3 s.d. 



CLs+b corresponds to the p- value for the s + b hypoth- 
esis. Figure 8 shows the quantity 1 — CLt, which is the 
p-value for the b hypothesis. These probabilities are local 
p-values, corresponding to searches for each value of Mh 
separately. These two p-values {CLg+b and 1 — CLh) 
provide information about the consistency of their re- 
spective hypotheses with the observed data at each value 
of Mh- Small values indicate rejection of the hypoth- 
esis and values above 50% indicate general agreement 
between the hypothesis in question and the data. As 
can be seen in Fig. 7, the observed value of CLg+b drops 
to sa 1% for Mh = 160 GeV, indicating limited consis- 
tency with the s + b hypothesis around this mass. In 
contrast, the observed value of CLs+b is close to unity 
for 120 < Mh < 145 GeV, whereas 1 - CU is small. 
At Mh = 125 (140) GeV, the value of 1 - CU is 4.1% 
(1.8%), corresponding to 1.7 (2.1) s.d. above the back- 
ground prediction. 

As a further investigation of this excess, we present 
in Fig. 9 the best fit of the data to the ratio of an to 
the SM prediction (crFit/csM)- The result of this fit, 
shown along with its band of ±1 s.d., yields a signal 
rate of approximately a factor of 1.4 larger than the SM 
cross section for Mh between 120 GeV and 145 GeV. For 
Mh = 125 GeV, we obtain a ratio of 1.4_± 0.9. The as- 
sociated production analyses with H ^ bb decay and the 
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FIG. 8: (color online) The observed (black solid line) and ex- 
pected (black short-dashed line) 1 — CLf, (background p-value) 
for all analyses combined for the range 90 < Mh < 200 GeV. 
Also shown is the expected background p-value for a presence 
of a Mh ~ 125 GeV SM Higgs boson signal in the data (blue 
long-dashed line). The shaded bands correspond, respectively, 
to the regions enclosing ±1 and ±2 s.d. fluctuations of the 
background. The three red horizontal dashed lines indicate 
the p-values corresponding to significances of 1, 2 and 3 s.d. 
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D0, L. , < 9.7 fb"^ 
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MH(GeV) 

FIG. 9: (color online) The best fit of the ratio (Jh /{o-h)bm 
as a function of Mh for all analyses combined for the 110 < 
Mh < 200 GeV. This indicates the values of the Higgs boson 
cross section that best match the observed data. The light 
shaded band indicates the ±1 s.d. region departure from the 
fit. The fit result is zero for masses below 110 GeV. 



H — !• W~^W~ — >■ (.'^vi~T' analyses dominate our sensitiv- 
ity. The dijet invariant mass resolution is approximately 
15% for associated production with H ^ hh decay. The 
mass resolution for the analyses with H — W^W~ de- 
cay is poor due to the undetected neutrinos in the final 
state. We therefore expect a Higgs boson signal to ap- 
pear as a broad excess over background, rather than a 
narrow resonance such as that expected at the LHC in 
the H ^ or H ^ ZZ 4£ final states. 

We study the excess at low mass by separating the 
major contributing sources according to the Higgs bo- 
son decay: H ^ bb, H ^ W^W' , H t+t^ and 
if — >■ 77 final states. Figure 10 shows the LLR values 
from the combination of the results from the ZH — )- Ubb, 
ZH — >■ uubb and WH — > tvbb searches, and illustrates a 
small excess that is compatible with the SM Higgs bo- 
son expected rate for 120 < Mh < 145 GeV. Figure 11 
shows the LLR values from the combination of the re- 
sults from searches for H — )• W'^W' — £v^v, H — 
W+W^ ivjj, and VH VWW, together with the 
VFVF-dominated subchannels from the H + X ^ iVijj 
analysis, and shows a similar excess of data over the back- 
ground for 110 < Mh < 150 GeV. At higher masses, 
where the Tevatron sensitivity to Higgs boson produc- 
tion is the largest, the LLR favors the b hypothesis. Fig- 
ure 12 shows the LLR values from the combination of 
the TT-dominated H + X —i- ir^jj subchannels and the 
VH — > ThThfJ. + X analysis, in which a significant frac- 
tion of the Higgs boson decays are to r+r^ pairs. Fig- 
ures 13-15, as well as Tables IV-VI, show the expected 
and observed 95% C.L. cross section limits in terms of 
ratio to the SM predictions 

for H ^ bb, H ^ W+W-, and H -> t+t" final 
states, respectively. The corresponding figures for the 
H ^ -f^ analysis can be found in Rcf. [30]. Figure 16 
shows the best fit of the ratio an ■ B/{(Th ■ B)sm for 



Mh = 125 GeV in each of the Higgs boson decay channels 
considered, as well as the central value for all analyses 
combined. These values are also given in Table VII. 
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FIG. 10: (color online) The observed (black solid line) and 
expected LLRs for the b (black short-dashed line) and s + b 
hypotheses (rod/light short-dashed line) for the combined 
WH/ZH,H bb analyses for the range 90 < Mh < 
150 GeV. The shaded bands correspond, respectively, to the 
regions enclosing ±1 and ±2 s.d. fluctuations of the back- 
ground. 
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FIG. 11: (color online) The observed (black solid line) and 
expected LLRs for the b (black short-dashed line) and s + b 
hypotheses (red/light short-dashed line) for the combined 
WH/ZH/H,H W+W^ analyses for the range 100 < 
Mh < 200 GeV. The shaded bands correspond, respectively, 
to the regions enclosing ±1 and ±2 s.d. fluctuations of the 
background. 



Interpretation in fourth generation and 
Fermiophobic Higgs boson models 



We also interpret our Higgs boson searches in models 
containing a fourth generation of fcrmions, and models 
with a fermiophobic Higgs boson. The fourth generation 
models [87] feature a modified Hgg coupling, leading to 
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TABLE IV: Expected (median) and observed upper limits for cr x B(H — bb) relative to the SM at 95% C.L. for the combined 
WH/ZH, H ^bb analyses for the range 90 < Mr < 150 GeV. 



Mh (GeV) 


90 


95 


100 


105 


110 


115 


120 


125 


130 


135 


140 


145 


150 


Expected 


1.29 


1.40 


1.21 


1.31 


1.45 


1.63 


1.92 


2.33 


2.99 


3.96 


5.52 


7.91 


11.35 


Observed 


0.96 


0.89 


1.05 


1.33 


1.51 


2.25 


2.96 


3.49 


4.29 


6.92 


8.65 


13.85 


13.90 



TABLE V: Expected (median) and observed upper limits for a x B{H W+W~) relative to the SM at 95% C.L. for the 
combined WH/ZH/H,H ^ W+W" analyses for the range 100 < Mr < 200 GeV. 



Mh (GeV) 
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115 
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7.25 
9.98 


7.09 
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6.49 
12.38 


5.34 
7.70 


3.97 
5.84 


2.92 
4.55 


2.33 
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1.88 
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1.64 
3.14 


1.40 
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1.20 
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1.02 
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0.75 
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0.66 
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FIG. 12: (color online) The observed (black solid line) and 
expected LLRs for the b (black short-dashed line) and s + b 
hypotheses (red/light short-dashed line) for the combined 
VH — >■ ThTh^J. + X and H -\- X ^ f-Thjj analyses for the range 
100 < Mh < 150 GeV. The shaded bands correspond, respec- 
tively, to the regions enclosing ±1 and ±2 s.d. fluctuations of 
the background. 



FIG. 13: (color online) Expected (median) and observed 
ratios for the upper limits of the cross section an at 
95% C.L. relative to the SM values for the combined 
WH/ZH,H bb analyses for the range 90 < Mh < 
150 GeV. The shaded bands correspond to the regions en- 
closing ±1 and ±2 s.d. fluctuations of the background, re- 
spectively. 



a nearly order of magnitude enhancement in the GGF 
cross section relative to the SM [88-90]. Previous inter- 
pretations of SM Higgs boson searches within the con- 
text of a fourth generation of fermions at the Fermilab 
Tevatron CoUider exclude 131 < Mr < 207 GeV [11]. 
Both ATLAS [91] and CMS [92] have performed sim- 
ilar searches, which exclude, respectively, 140 < Mh < 
185 GeV and 110 < Mh < 600 GeV. Ahhough the larger 
coupling increases the decay width to gg, the WW* de- 
cay mode remains dominant for Mh > 135 GeV. There 
is also a small contribution from H ZZ* — Uvi/ 
production that increases with Mh- Wc consider two 
fourth generation scenarios: (i) a "low mass" scenario in 
which the mass of the fourth generation neutrino is set 
to m^4 = 80 GeV, and the mass of the fourth generation 
charged lepton m£4 is set to 100 GeV, and (ii) a "high 
mass" scenario in which m^4 = mn = 1 TeV, so that the 
fourth generation leptons do not affect the decay branch- 
ing fractions of the Higgs boson. In both scenarios the 



fourth generation quark masses are set to be those of the 
high mass scenario in Ref. [90] . 

Wc consider only gg ^ H production and the H — > 
W+W- l+iyl-D and H -> W+W~ tvq'qim) chan- 
nels to set limits on the fourth generation models, and 
also set a limit on a{gg H) x B{H W+W-). We 
scale the product of the cross sections and branching frac- 
tions to the results from hdecay, modified to include 
the fourth generation. We retrain our multivariate dis- 
criminants to take only the above signals into account, 
and do not include events with two or more jets in the 
H — > W^W~ — >■ evev / fsv fiv analyses. We also do not in- 
clude the theoretical uncertainty on (j{gg — >■ H) x B{H — > 
W+iy~) since the absolute cross section limits do not 
depend on the prediction. We include the theoretical un- 
certainties for limits on ratios to cross sections. 

Figure 17 shows the combined limits on (j{gg — > 
H) X B{H W^W~), along with the fourth genera- 
tion theory predictions for the high mass and low mass 
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TABLE VI: Expected (median) and observed upper limits for a x B{H — >• r+r ) relative to the SM at 95% C.L. for the 
combined VH — >■ ThThfi + X and + X — >■ irhjj analyses for the range 100 < Mh < 150 GcV. 



Mh (GeV) 
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150 


Expected 
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6.54 


6.59 
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8.46 


9.05 
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11.28 


12.11 


Observed 


8.42 
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18.59 
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FIG. 14: (color online) Expected (median) and observed 
ratios for the upper limits of the cross section an at 
95% C.L. relative to the SM values for the combined 
WH/ZH/H,H W+W analyses for the range 100 < 
Mh < 200 GeV. The shaded bands correspond to the regions 
enclosing ±1 and ±2 s.d. fluctuations of the background, re- 
spectively. 
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FIG. 16: (color online) The best fit of an ■ B/{aH ■ B)sm for 
various Higgs boson decays for Mh = 125 GeV. The central 
value for all combined analyses is shown with its 1 s.d. band 
(shaded area). 
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TABLE VII: The best fit Higgs boson cross section times 
branching fraction as a ratio to the SM cross section times 
branching fraction for AIh = 125 GeV for the individual com- 
binations according to Higgs boson decay mode, as well as the 
full combination. 



Combined 



1.40 



+0.92 

0.88 



H 
H 
H 
H 



w+w- 

■T+T- 

hb 



gQ + 1.63 

3.96 
1.23 



1.52 

+4.11 
3.38 

+ 1.24 



FIG. 15: (color online) Expected (median) and observed 
ratios for the upper limits of the cross section an at 
95% C.L. relative to the SM values for the combined VH — >■ 
ThThf^ + X and H + X ^ t-Thjj analyses for the range 
100 < Mh < 150 GeV. The shaded bands correspond to 
the regions enclosing ±1 and ±2 s.d. fluctuations of the back- 
ground, respectively. 



scenarios. Wc exclude a SM-like Higgs boson in the 
range 125 < Mr < 218 GeV at 95% C.L., with an ex- 
pected exclusion range of 122 < Mh < 232 GeV in the 
low mass scenario. In the high mass scenario, the ob- 
served (expected) exclusion range is 125 < Mh < 228 
(122 < Mh < 251) GcV. 

In the fermiophobic model (FHM), the lightest Higgs 
boson Hf couplings to fermions vanish at leading order, 
but otherwise Hf is like the SM Higgs boson. Hence, 
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gg — )• Hf production is negligible, and Hf decays to 
fermions are forbidden, but V + Hf and vector boson 
fusion qq — > q'q'Hf production remain nearly unchanged 
relative to the SM. The WW, ZZ, 77, and Zj decays 
comprise nearly the entire decay width. For all the 
Hf — > W^W" decay has the largest branching fraction. 
The Hf — >■ 77 branching fraction is greatly enhanced 
over the SM for all Mhj, and it provides most of the 
search sensitivity for Mh^ < 120 GeV. 

The CDF and DO Collaborations have previously pub- 
lished results in the Hf — > 77 decay channel [93, 94]. The 
analyses described here supersede previous FHM searches 
at DO. The ATLAS and CMS CoUaborations have 
performed fermiophobic searches, and exclude 110 < 
MHf < 118.0 GeV, 119.5 < A/^^ < 121.0 GeV [95], 
and 110 < MHf < 147 GeV [92] using 77 final states, 
and 110 < < 194 GeV when other final states are 

included [96]. 

We combine the 77 and H W^W~ decay 

channels, produced either in association with a Y boson, 
or in VBF, for the FHM interpretation. We reoptimize 
the SM 77 — >■ 77 analysis to take into account the dif- 
ferent kinematics in the FHM, e.g., the presence of an 
associated vector boson in the FHM, or recoiling quark 
jets in VBF, which shift the transverse momentum spec- 
trum of the Higgs boson to higher values than in the 
SM. Likewise, we retrain the multivariate discriminants 
for the H W^W~ —J- (e+e", /z"''^", e^^^)i^i? anal- 
yses to account for the suppressed GGF process in the 
FHM. We retain the existing subdivision into categories 
that are based on the number of reconstructed jets in the 
event. The other SM H — >■ WW analyses can be inter- 
preted directly in the FHM without reoptimization, af- 
ter separating the relative contributions from GGF, WH, 
ZH, and VBF in each contributing channel, removing the 
GGF component, and scaling the remaining signal con- 
tributions by the ratio of the branching fraction in the 
FHM and SM, B{Hf ^ WW)/B{Hsm ^ WW). Fig- 
ure 18 shows the combined FHM limits. The observed 
(expected) 95% C.L. exclusion range is 100 < MHf < 114 
(100 < MHf < 117) GeV. 



VIII. CONCLUSIONS 
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FIG. 17: (color online) Expected and observed 95% C.L. up- 
per limits on Higgs boson production in fourth generation 
models as a function of Higgs boson mass. The blue and red 
lines represent the theoretical predictions with its uncertain- 
ties in the fourth generation "high mass" and "low mass", 
respectively. Below 160 GeV the models overlap and have 
similar uncertainties. When setting these limits we do not in- 
clude the theoretical cross section uncertainties. The shaded 
bands correspond to the regions enclosing ±1 and ±2 s.d. fluc- 
tuations of the background. 
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FIG. 18: (color online) Expected and observed 95% C.L. up- 
per limits on fermiophobic Higgs boson production as function 
of Higgs boson mass. We exclude a fermiophobic Higgs boson 
with a mass below 114 GeV. The shaded bands correspond 
to the regions enclosing ±1 and ±2 s.d. fluctuations of the 
background. 



We have presented a combination of searches for SM 
Higgs boson production with the DO experiment using 
data corresponding to up to 9.7 fb~^ of pp collisions at 
— 1.96 TcV. Wc set upper limits on the production 
cross section at 95% C.L. for Higgs boson masses of 90 < 
Mh < 200 GeV. We also interpret the searches in terms 
of models containing a fourth generation of fermions, as 
well as models with a fermiophobic Higgs boson {Hf) 
having suppressed couplings to fermions. We exclude 
a Higgs boson in the mass range 125 < Mh < 218 
(125 < Mh < 228) GeV, in the low mass (high mass) 
fourth generation scenario, and a fermiophobic Higgs 



boson with a mass 100 < MHf < 114 GeV. The ob- 
served upper limits on SM Higgs boson production are 
2.86 (0.66) X asM at Mh = 125 (165) GeV, with an ex- 
pected limit of 1.68 (0.70) x asm- We exclude the regions 
of 90 < Mh < 101 GeV and 157 < Mh < 178 GeV with 
an a priori expected exclusion of 155 < Mh < 175 GeV. 
In the range of Mh ~ 120—145 GeV, the data exhibit an 
excess above the background prediction of up to two stan- 
dard deviations consistent with the presence of a 125 GeV 
SM Higgs boson. Each of the four main Higgs boson de- 
cay mode combinations contributes to this excess. The 
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analyses combined here also provide inputs to the overall 
Tevatron combination, which shows evidence for a SM- 
like Higgs boson [97]. 
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